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PREFACE 

The  study  presented  herein  was  conducted  by  the  Arnold  Engineering  Development 
Center  (AEDC),  Air  Force  Systems  Command  (AFSC),  under  Program  Element  63723F. 
This  work  was  done  by  ARO,  Inc.,  (a  subsidiary  of  Sverdrup  &  Parcel  and  Associates, 
Inc.),  contract  operator  of  AEDC,  AFSC,  Arnold  Air  Force  Station,  Tennessee.  The 
research,  conducted  under  ARO  Project  Number  V37A-32A,  was  done  in  support  of  the 
High  Reynolds  Number  Wind  Tunnel  (HIRT)  Project.  The  author  of  this  report  was 
Frederick  L.  Shope,  ARO,  Inc.  The  manuscript  (ARO  Control  No.  ARO-VKF-TR-75-72) 
was  submitted  for  publication  on  March  17,  1975. 

The  author  wishes  to  acknowledge  the  contributions  of  Rogers  F.  Starr  who  carefully 
supervised  and  steered  the  effort  and  made  many  contributions  to  the  experiment  as  well 
as  the  theory. 
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1.0  INTRODUCTION 

This  study  was  part  of  the  effort  to  develop  a  transonic,  high  Reynolds  number 
wind  tunnel  (HIRT)  based  on  the  Ludwieg  tube  concept.  The  proposed  HIRT  Facility 
(Ref,  1)  consisted  of  a  long  slender  pressure  bottle  (15-ft-diameter  by  1 ,660-ft-long)  which 
would  provide  high-pressure  (500-psi  stagnation)  and  low-temperature  (-60°F  stagnation) 
air  to  a  porous  walled,  transonic  test  section.  To  achieve  these  stagnation  conditions  in 
the  test  section  during  a  tunnel  run,  the  air  in  the  supply  tube  would  have  to  be  compressed 
to  700  psi  and  cooled  to  -30°F.  In  order  for  the  charge  air  to  remain  at  the  desired 
temperature  until  a  tunnel  run  could  be  made,  the  steel  supply  tube  itself  would  have 
to  be  cooled  to  -30“F.  In  addition,  the  cooling  would  have  to  be  accomplished  in  a 
reasonable  length  of  time,  and  the  final  temperature  distribution  would  have  to  be  uniform 
(within  5°F).  One  proposed  system  to  achieve  these  conditions  would  consist  of  a  series 
of  compressors  and  refrigeration  units  to  recirculate  very  cold  air  (-60°F  to  -80°F)  through 
the  supply  tube  to  cool  it  by  forced  convection.  However,  the  flow  rates  (<200  Ib/sec) 
practically  attainable  from  reasonably  sized  compressors  and  the  high  pressures  (~700  psi) 
required  to  avoid  inordinately  large  recirculation  and  refrigeration  plumbing  would  also 
result  in  very  low  bulk  tube  velocities  and  large  densities.  Hence,  the  flow  could  well 
be  dominated  by  free  convection  and  induce  unacceptably  large  temperature  gradients 
in  the  tube.  In  addition,  the  size  of  the  compressors  and  refrigeration  units  would  have 
a  first  order  impact  upon  the  total  facility  cost.  Therefore,  the  current  work  was  undertaken 
to  determine  if  the  supply  tube  could  be  suitably  cooled  by  simply  circulating  cold  air 
from  one  end  to  the  other. 

This  report  documents  the  results  of  a  study  of  internal  convection  cooling  of  a 
long  slender  tube  in  which  the  classical  free  and  forced  convection  heat-transfer  mechanisms 
were  severely  disrupted  by  a  high-speed  inlet  jet.  The  flow  regime  of  the  test  corresponded 
to  a  tube  diameter  Reynolds  number  of  4  x  10“^  and  Grashof  numbers  varying  from  10^ 
to  10^0  (Grp /Rep  ~  10  to  100).  These  values  resulted  from  relatively  low  velocities, 
approximately  0.2  to  0.7  ft/sec,  and  high  pressures  6  to  28  atmospheres. 

The  experiment  was  performed  using  a  steel  tube  about  1-in.  thick.  111  feet  long, 
and  14-in.  ID.  Exact  dimensions  are  listed  in  Table  1.  Air,  cooled  to  between  -5°F  and  20°F 
by  Joule-Thompson  expansion  from  a  4,000-psi  air  supply,  was  injected  into  one 
end  of  the  tube  through  a  relatively  high-speed  jet  (200  to  300  ft/sec)  normal  to  the 
axis  of  the  tube.  Data  collected  included  the  axial  and  circumferential  timewise  temperature 
distributions  in  the  steel  tube  and  the  axial  and  vertical  timewise  temperature  distributions 
of  the  air. 

Two  things  make  the  current  experiment  somewhat  unique  among  heat-transfer  tests. 
First,  the  relatively  high  pressures  and  low  bulk  velocities  place  the  physical  situation 
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at  the  extremes  of  the  Reynolds  number  and  Grashof  number  regimes  usually  covered 
by  the  classical  heat-transfer  experiments.  However,  the  low  Reynolds  numbers  result  not 
from  high  viscosity,  as  in  the  oil-flow  experiments,  but  from  low  velocities,  and  the  high 
Grashof  numbers  from  high  pressure  rather  than  size  of  the  apparatus.  A  second  unique 
characteristic  of  this  experiment  is  that  the  tube  was  insulated  from  the  environment  and 
its  temperature  was  allowed  to  change  locally  in  accordance  with  the  heat  transfer  required 
by  the  fluid  flow.  This  is  in  contrast  to  the  classical  experiments  where  the  tube  is  usually 
provided  an  external  energy  source  to  keep  it  at  a  constant  temperature.  Allowing  the 
tube  temperature  to  vary  temporally  and  spatially  is,  of  course,  necessary  in  studying 
the  development  of  temperature  gradients. 

The  report  documents  the  experimental  apparatus,  test  procedure,  data  reduction  and 
interpretation,  analysis,  and  conclusions.  Details  of  the  data  reduction  and  modeling  are 
given  in  the  Appendixes. 


2.0  APPARATUS 

The  experimental  apparatus  was  adapted  from  the  supply  tube  of  the  pilot  HIRT 
in  the  von  Karman  Gas  Dynamics  Facility  (VKF)  (Ref.  1).  The  characteristics  of  the  tube 
are  listed  in  Table  1.  To  prepare  the  tube  for  the  test,  it  was  disconnected  from  the 
test  section,  though  a  short  (18.5-in.)  piece  of  the  contraction  section  (a  converging  nozzle) 
was  left  attached  to  the  tube  in  order  to  use  an  existing  blind  flange  to  seal  the  pressure 
vessel.  To  measure  the  temperature-time  histories  and  spatial  distributions  in  the  tube  steel, 
the  tube  was  initially  instrumented  with  18  thermocouples  as  illustrated  in  Fig.  1  (open 
squares).  At  each  of  three  stations  along  the  tube,  seven  thermocouples  were  installed 
as  shown  in  Section  AA.  Five  copper-constantan  thermocouples  were  spaced  uniformly 
on  the  outside  of  the  tube,  and  one  was  potted  into  a  hole  in  the  tube  to  measure  inside 
surface  temperature.  In  addition,  at  each  of  these  three  stations,  a  Chromel®-constantan 
thermocouple  on  the  end  of  a  vertically  sliding  tube  was  installed  to  measure  vertical 
gradients  in  the  air. 


Table  1.  Tube  Characteristics 


Length 


111  ft 


Inside  diameter 


13.94  in. 


Thickness 
Density 
Specific  heat 
Thermal  conductivity 


1.026  in. 

0.283  Ibm/in.^ 
0.11  Btu/lbm 
27  Btu/hr-ft-OR 


Inside  roughness  0.0004-in.  to  0.0007-in.  roughness 

superimposed  on  a  larger  scale  0.001-in. 
to  0.0015-in.  roughness  of  more  or  less 
flat  plateaus  and  valleys  (-0.25  in. 
across) 
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o  ♦  Thermocouples  Measuring  Air  Temperature 
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Section  •7.34 
in.  by  9. 17  in. 
Rectangular 
Exiti 


Figure  1.  Schematic  of  experimental  apparatus. 


The  system  to  handle  the  coohng  air  from  the  4,000-psi  bottle  included  an  automatic 
control  valve  which  throttled  the  air  to  a  low  enough  pressure  to  prevent  choking  at 
the  inlet  port.  Prior  to  entering  the  tube,  the  3/4-in.  supply  hne  was  expanded  into  a 
2.5-in.-diam  stilling  chamber  where  the  stagnation  temperature  of  the  inlet  air  was 
measured.  The  air  then  passed  directly  from  the  stilling  chamber  into  the  tube  through 
a  few  inches  of  3/4-in.  pipe.  Air  was  exhausted  from  the  downstream  end  of  the  tube 
through  a  control  valve  and  sonic  orifice  flowmeter  (Ref.  2),  an  arrangement  used  to 
control  and  measure  the  mass  flow  rate  through  the  tube.  The  exhaust  line  leaving  the 
tube  was  instrumented  with  a  thermocouple  to  measure  the  exit  air  temperature  upstream 
of  the  control  valve.  The  flowmeter  also  contained  a  thermocouple  upstream  of  the  orifice 
as  well  as  pressure  gages  on  either  side  of  the  orifice  plate  to  ensure  the  pressure  drop 
across  the  plate  was  sufficient  to  keep  it  choked.  Downstream  of  the  meter  the  air  was 
exhausted  to  atmosphere.  Pressure  in  the  tube  was  measured  with  a  precision  pressure 
transducer  located  approximately  at  the  midpoint  of  the  tube.  The  very  low  velocities 
and  correspondingly  small  pressure  drops  along  the  tube  did  not  warrant  the  installation 
of  additional  pressure  taps  along  the  tube.  Since  calculations  indicated  that  free  convective 
heat  transfer  from  the  outside  of  the  tube  was  important,  the  entire  tube  was  wrapped 
in  fiber  glass  insulation  backed  with  aluminized  Mylar®  to  reduce  losses.  The  inlet  and 
outlet  plumbing  were  wrapped  with  insulation  tape. 


Three  cooling  runs  were  made  with  this  apparatus,  and  the  results  indicated  the  need 
for  more  detailed  information  on  the  development  of  circumferential  and  axial  temperature 


7 


AEDC-TR-75-142 


gradients  in  the  tube  as  well  as  the  possibility  of  reducing  these  gradients  with  additional 
inlet  ports.  To  collect  this  additional  temperature  gradient  data,  the  tube  was  further 
instrumented  with  thermocouples  at  eleven  additional  stations  along  the  tube  (black  squares 
in  Fig.  1).  At  each  station  a  thermocouple  was  attached  at  the  top  and  bottom  of  the 
tube.  In  addition,  a  second  inlet  port  and  stilling  chamber  were  installed  about  15  diameters 
downstream  of  the  first  inlet  port.  The  plumbing,  valves,  and  pressure  gages  necessary 
for  control  of  the  two  inlet  ports  were  also  added  to  the  system. 

The  data  collection  was  performed  with  a  slow  speed  scanning  digital  voltmeter  which 
scanned  the  25  to  50  channels  of  thermocouples  and  pressure  taps.  The  first-line  data 
reduction  was  performed  during  the  test  using  a  digital  computer  to  convert  the  voltage 
printouts  to  temperatures. 


3.0  PROCEDURE 

Six  separate  cooling  tests  were  begun  by  first  pumping  the  charge  tube  up  to  the 
desired  pressure  with  the  exit  valve  closed.  Upon  reaching  the  charge  pressure,  the  inlet 
and  outlet  valves  were  adjusted  to  produce  the  desired  mass  flow  rate  and  hold  the  charge 
pressure  constant.  At  intervals  from  15  to  30  minutes,  readings  were  taken  from  all  data 
channels.  Measurements  were  taken  from  the  three  air  survey  probes  at  five  vertical 
positions  by  manually  sliding  the  probes  up  and  down.  At  the  conclusion  of  the  test, 
the  flanges  capping  the  ends  of  the  tube  were  removed  to  allow  the  tube  to  warm  up 
to  room  temperature  for  the  next  test. 

On  certain  individual  tests,  various  special  operations  were  performed.  On  Run  No. 
1 ,  following  completion  of  the  cooling  process,  the  tube  was  allowed  to  stand  idle  at 
pressure  with  no  mass  flow  to  determine  how  rapidly  the  circumferential  and  axial 
temperature  gradients  in  the  tube  would  dissipate.  On  Run  No.  2,  this  was  tried  at 
atmospheric  pressure  rather  than  charge  pressure.  On  Run  No.  4,  the  effects  of  staging 
on  the  cooling  process  were  measured.  On  this  run,  cooling  air  was  injected  through  the 
first  air  inlet  port  for  several  hours  until  the  steel  in  the  region  of  the  inlet  approached 
the  inlet  air  temperature.  The  first  inlet  port  was  then  closed  and  the  second  one  opened 
(15  diameters  downstream)  for  the  remainder  of  the  test. 

Sonic  orifice  calibrations  were  performed  by  pumping  the  tube  to  a  given  pressure, 
shutting  off  the  inlet  air,  and  then  dumping  the  tube  through  the  sonic  orifice  while 
measuring  pressure  and  temperature-time  histories  in  the  tube  and  sonic  orifice. 
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4.0  RESULTS 

The  data  to  be  discussed  below  were  obtained  after  two  levels  of  data  reduction. 
The  first  level,  performed  during  the  test,  consisted  primarily  of  converting  the 
thermocouple  voltages  to  temperatures  using  either  curve-fits  of  Ref.  3  or  curve-fits  to 
the  standard  tabular  data.  These  temperature-time  histories  were  later  converted  to 
heat-transfer  coefficients  and  Nusselt  numbers  according  to  the  computational  procedure 
discussed  in  Appendix  A.  The  conditions  for  each  of  the  six  runs  are  summarized  in  Table 
2. 


Table  2.  Test  Conditions 


Hun  Number 

- 1 

2 

3 

4 

5 

6 

Charge  Pressure^  psig 

400 

150 

75 

1  50 

75 

150 

Average  Mass  Flo«  Rate,  Ibm/sec 

0.41 

0.37 

0 . 34 

0.33 

0.33 

0.33 

Test  Time,  Air  Floning,  hr 

4 .9 

7  .0 

5.0 

7 . 5 

3.9 

3  .0 

Idle  Period  at  End  of  Test,  hr 

0.9  at  pressure 

1.0  at  ambien  t 

0 

0 

0 

0 

Inlet  Temperature  Range, 

475  to  490 

454  to  475 

462  to  482 

450  to  477 

460  to  481 

460  to  488 

Average  Velocity  In  Charge  Tube,  ft/sec 

0.  17 

0 . 40 

0.66 

0  34 

0.65 

0 . 36 

Reynolds  Number  Based  on  Tube  Diameter 

39,000 

36,000 

33,000 

33,000 

32 , 000 

31 . 000 

Grashof  Number  Based  on  Tube  Diameter 

2  X  lolO 

5  X  10* 

2  X  10® 

5  X  10® 

2  X  10® 

5  X  10® 

Initial  Tube  Temperature  Range,  ‘^R 

525  to  529 

524  to  531 

526  to  529 

530  to  536 

528  to  530 

530  to  534 

Final  Tube  Temperature  Range, 

496  to  SIS 

475  to  503 

486  to  517 

477  to  507 

479  to  522 

468  to  522 

Prandtl  Number 

0.68 

0.68 

0.68 

0  69 

0.6B 

0.68 

Some  of  the  temperature-time  histories  from  Runs  No.  1,  2,  and  3  are  presented 
in  Figs.  2  through  6.  For  these  first  tests,  temperature  measurements  were  limited  to 
three  axial  stations  along  the  length  of  the  tube  (Fig.  1)  located  roughly  near  the  inlet, 
the  midpoint,  and  the  outlet.  Tests  were  conducted  at  pressures  of  400,  150,  75  psig 
for  nominally  constant  mass  flow  rates  of  0,3  to  0.4  Ib/sec.  Figure  2  shows  how  the 
steel  temperature  varied  during  six  hours  of  Run  No.  2  (150  psig).  As  expected,  the  steel 
temperature  dropped  most  rapidly  near  the  inlet  since  this  is  the  region  of  greater  air-steel 
temperature  difference  early  in  the  test.  The  two  stations  further  downstream  showed 
correspondingly  lesser  slopes  since  they  were  cooled  by  increasingly  warmer  air.  For  each 
station,  curves  are  shown  for  both  the  top  and  bottom  outside  temperature.  At  Station 
1,  there  was  essentially  no  difference  in  temperature  between  the  top  and  bottom  of 
the  tube.  At  Station  2  there  was  a  marked  difference  -  as  great  as  10°R  at  3.5  hours 
-  and  Station  3  showed  a  significant  but  lesser  difference.  The  explanation  of  these 
experimental  results  is  believed  to  lie  in  accounting  for  the  turbulence  caused  by  the  inlet 
jet,  which  reached  velocities  of  200  to  300  ft/sec.  After  presentation  of  the  remaining 
experimental  results,  the  turbulent  mechanism  believed  responsible  will  be  discussed  in  more 
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detail  along  with  an  attempt  to  model  these  results.  One  additional  curve  on  Fig.  2  shows 
how  the  inlet  air  temperature  varied  during  the  test.  The  20‘‘R  variation  shown  here  is 
attributable  to  variation  of  the  pressure  of  the  air  supply,  which  was  nominally  3,500 
psig  but  occassionally  dropped  below  2,000  psig.  The  variation  of  the  supply  pressure 
was  due  to  other  heavy  users  of  the  supply  bottle  and  resulted  in  varying  Joule-Thorn pson 
drops  in  stagnation  temperature.  The  effect  of  the  inlet  temperature  rise  beginning  at 
about  one  hour  can  be  seen  in  the  slightly  lower  slope  of  the  steel  temperature  between 
two  and  three  hours,  and  the  gradual  warming  at  Station  1  beyond  4.5  hours  resulted 
when  the  inlet  air  temperature  actually  rose  above  the  local  steel  temperature.  With  the 
exception  of  the  top-to-bottom  temperature  differences,  the  data  shown  in  Fig.  2  for 
Run  No.  2  may  be  considered  typical  of  all  six  runs. 


Time,  hr 


Figure  2.  Steel  temperature  as  a  function  of  time 
for  Run  No.  2  (150  psig). 

The  circumferential  temperature  difference  in  the  tube  is  of  course  due  to  vertical 
temperature  gradients  in  the  cooling  air.  Figures  3,  5,  and  6  show  the  data  obtained  with 
the  sliding  air  survey  probes.  Figure  3  compares  the  vertical  gradients  at  each  axial  station 
for  Run  No.  1  (400  psig).  Note  that  significant  gradients  occur  at  all  three  stations,  but 
that  the  gradients  tend  to  decrease  in  the  downstream  direction.  Figure  4  shows  similar 
data  for  the  circumferential  steel  temperature.  Note  the  large  gradients  at  Stations  1  and 
2  during  most  of  the  run.  The  gradients  at  Station  3  before  the  two-hour  point 
are  small  but  consider  also  that  little  cooling  has  taken  place  (the  air  is  relatively  warm 
by  the  time  it  reaches  Station  3).  As  Station  3  begins  to  cool  beyond  two  hours,  however. 
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gradients  do  begin  to  appear.  Figure  5  shows  data  for  150  psig,  and  there  are  marked 
differences  from  the  400-psig  case.  There  are  essentially  no  gradients  at  Station  1  while 
Station  2  shows  heavy  stratification  and  Station  3  somewhat  lesser  gradients  as  compared 
to  Fig.  5b.  Figure  6  (75  psig)  shows  no  stratification  at  the  first  station,  slight  but 
temporary  stratification  at  the  second,  and  significant  stratification  at  the  third. 


Sym  Time,  hr 
"  0 

Axial  Slation  1 


Axial  Station  2 


Figure  4.  Circumferential  steel  temperature  distribution  for 
Run  No.  1  (400  psig). 


At  the  completion  of  Runs  No.  1  and  2,  after  the  cooling  air  had  been  shut  off, 
the  tube  was  allowed  to  stand  idle  for  an  hour  while  the  gradients  dissipated.  The  idle 
period  following  Run  No.  1  was  at  full-charge  pressure  of  400  psig,  while  Run  No.  2 
was  at  ambient  pressure.  Figure  7  shows  how  the  top-to-bottom  temperature  difference 
in  the  tube  varied  during  the  idle  hour  at  each  of  the  three  stations.  The  gradients  were, 
of  course,  much  larger  at  the  end  of  Run  No.  1  because  of  the  higher  operating  pressure; 
by  0.8  hours  after  shutdown,  the  gradients  had  dropped  to  less  than  5°R.  Following  Run 
No.  2  the  gradients  were  already  less  than  5°R  but  temporarily  increased  slightly  before 
dissipating,  due  probably  to  the  sudden  loss  of  the  jet  turbulence  and  setting  up  of  a 
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Strong  free  convection  cycle.  The  rise  in  gradients  did  not  appear  in  the  first  case  because 
they  were  already  so  large  and  the  jet  turbulence  was  of  little  importance  in  the 
high-pressure  case. 


Figure  7.  Dissipation  of  circumferential  temperature  gradients  in 
tube  following  shutdown  of  cooling  air. 


To  refine  the  results  of  the  first  three  runs  and  study  the  effects  of  stratification 
in  more  detail,  three  additional  runs  were  made  after  further  instrumenting  the  tube  as 
explained  previously.  Figures  8  and  9  show  how  the  top-to-bottom  temperature  gradients 


Figure  8.  Effect  of  stratification  for 
Run  No.  5  (75  psig). 


15 


AEDC-TR-75-142 


in  the  tube  varied  with  time  and  axial  position  along  the  tube.  At  75  psig  (Fig.  8)  the 
gradients  are  seen  to  be  essentially  zero  out  to  28  diameters  but  show  a  sudden  increase 
from  that  point  on.  At  twice  the  pressure  (Fig.  9),  the  small  gradients  extend  only  to 
24  diameters  and  then  increase  much  more  rapidly  from  there  on.  Run  No.  4  (150  psig) 
was  performed  to  measure  the  effect  upon  these  gradients  of  moving  the  inlet  jet  further 
downstream  once  an  initial  section  of  the  tube  had  been  cooled  significantly  below  ambient. 
At  5.5  hours  into  the  run,  air  inlet  port  No.  1  (Fig.  1)  was  closed  and  port  No.  2  was 
opened.  Figure  10  shows  the  large  axial  temperature  gradients  in  the  tube  just  before 


Axial  Station,  diameters 


Figure  9.  Effect  of  stratification  for  Run  No.  6  (150  psig). 


Figure  10.  Steel  temperature  at  top  of  tube  versus  axial  station  for 
Run  No.  4  (150  psig)  before  and  after  staging. 
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staging  (changing  of  inlet  ports)  at  5.5  hours.  Two  hours  later  the  gradients  upstream 
of  the  second  inlet  had  essentially  washed  out.  Figure  1 1  shows  a  similar  effect  on  the 
circumferential  temperature  gradients.  Note  again  the  sudden  increase  in  the  gradients  after 
a  certain  distance  from  the  inlet. 


Axial  Station  Number 

Figure  11.  Difference  between  top  and  bottom  steel  temperatures  for 
Run  No.  4  (150  psig). 

From  the  temperature-time  histories  of  the  air  and  steel,  heat-transfer  coefficients 
along  the  tube  may  be  computed  as  outhned  in  Appendix  A.  Figure  12  shows  the  results 
of  this  calculation,  which  accounted  for  the  vertical  and  circumferential  temperature 
gradients  with  an  averaging  technique.  Gradients  in  the  radial  direction  (normal  to  the 
tube  wall)  were  shown  to  be  nearly  zero  by  the  thermocouples  mounted  on  the  inside 
wall  and  were  therefore  ignored  in  the  data  reduction.  In  Fig.  12,  the  experimental 
heat-transfer  coefficients  are  shown  as  scatter  bands.  The  scatter  at  88  diameters,  for 
instance,  can  be  charged  to  the  very  small  air-steel  temperature  differences  which  result 
when  the  air  has  been  warmed  to  near  the  local  tube  temperature.  Thus,  when  the  small 
net  heat  transfer,  as  manifested  by  the  small  timewise  change  in  steel  temperature,  is 
divided  by  the  air-steel  temperature  difference,  which  is  also  small,  the  uncertainty  in 
the  calculated  heat-transfer  coefficient  becomes  somewhat  large.  The  experimental  data 
are  compared  with  three  classical  Nusselt  number  relations:  the  turbulent  flat  plate 
boundary  layer,  the  fully  developed  turbulent  pipe,  and  free  convection  in  an  enclosed 
space  (Ref.  4).  It  is  clear  that  none  of  the  classical  relations  correlate  well  with  the 
experimental  results,  though  far  downstream  free  convection  gives  the  closest  estimate. 
In  fact,  the  heat-transfer  mechanism  occurring  in  this  experiment  produced  heat-transfer 
coefficients  which  are  two  to  five  times  greater  than  anything  expected  from  classicial 
theory. 
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Figure  12.  Experimental  heat-transfer  coefficients  compared  with 
classical  Nusselt  number  relations. 


In  summary,  the  following  experimentally  derived  facts  may  be  compiled. 

1.  The  flow  in  the  tube  stratifies  into  layers  of  warmer  and  cooler  air,  the 
degree  of  stratification  depending  upon  the  charge  pressure  and  distance 
from  the  inlet. 

2.  At  the  lower  pressures  and  sufficiently  close  to  the  inlet  port,  the 
stratification  is  essentially  zero  but  seems  to  grow  suddenly  at  some  specific 
distance  from  the  inlet,  also  depending  upon  the  pressure. 

3.  The  mechanism  responsible  for  these  phenomena  produces  rates  of  heat 
transfer  much  greater  than  expected  from  any  classicial  theory. 

5.0  DISCUSSION 

The  three  experimental  conclusions  listed  in  the  preceding  section  appear  to  be 
logically  explainable  in  the  following  terms.  The  inlet  jet  (200  to  300  ft/sec)  sets  up 
a  relatively  intense  turbulence  in  the  upstream  portion  of  the  tube.  This  turbulence 
immediately  begins  to  decay  in  intensity  through  viscous  dissipation  and  is  simultaneously 
convected  downstream  by  the  bulk  velocity  required  by  the  net  mass  flow  rate  through 
the  tube.  Nearer  the  inlet,  the  turbulence  is  sufficiently  intense  to  prevent  development 
of  either  a  regular  free  convection  flow  or  a  boundary  layer.  Further  downstream,  the 
turbulence  eventually  decays  to  the  point  where  it  no  longer  has  sufficient  intensity  to 
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prevent  buoyant  forces  from  inducing  stratication  of  the  flow  into  warm  and  cold  layers. 
Simultaneously  with  the  decay  and  convection  of  the  turbulence,  an  increasing  amount 
of  air  becomes  heated  above  the  inlet  temperature,  gradually  strengthening  the  importance 
of  buoyancy  relative  to  the  turbulence  as  the  warmed  masses  of  air  are  convected 
downstream.  The  effect  of  the  turbulence  on  the  heat-transfer  process  is  to  prevent  the 
development  of  any  form  of  boundary  layer  and  to  cycle  much  colder  air  within  a  distance 
from  the  wall  much  less  than  that  of  a  boundary-layer  thickness.  As  the  turbulence  decays 
and  buoyancy  increases,  free  convection  gradually  takes  over  as  the  dominant  heat -transfer 
process  with  stratification  as  the  consequence. 

In  an  attempt  to  correlate  this  proposed  flow  model  with  the  experimental  results, 
a  simplified  mathematical  model  has  been  developed  in  Appendix  B.  The  model  treats 
the  heat-transfer  process  as  one  of  pure  conduction  across  a  thin  layer  of  fluid  with  a 
thickness  of  the  order  of  a  laminar  sublayer  reduced  in  thickness  by  the  jet  turbulence. 
Comparison  of  the  laminar  sublayer  thickness  for  a  turbulent  pipe  flow  not  disrupted 
by  a  Jet  with  a  reduced  sublayer  needed  to  produce  the  experimental  heat-transfer 
coefficients  allows  an  approximate  estimate  of  the  turbulent  intensity  and  its  variation 
with  distance  from  the  inlet.  A  second  part  of  the  model  postulates  that  the  onset  of 
stratification  may  be  correlated  with  axial  variation  of  the  ratio  of  the  turbulent  energy 
to  the  buoyant  energy  and  that  when  this  energy  ratio  falls  below  a  critical  value, 
stratification  begins.  The  resulting  quantitative  model  predicts  that  the  energy  ratio  at 
any  given  axial  station  is  inversely  proportional  to  the  fourth  root  of  the  pressure,  which 
agrees  roughly  with  the  points  (on  Figs.  8  and  9)  at  which  stratification  occurs. 

The  model  described  above  is  not  the  only  possible  explanation  of  the  experimental 
results,  and  two  alternate  hypotheses  were  considered.  The  first  of  the  two  ideas  which 
received  considerable  attention  was  an  explanation  of  the  onset  of  stratification  in  terms 
of  a  growing  thermal  boundary  layer.  The  essential  idea  was  that  a  warm  boundary  layer 
grew  on  the  wall  and  merged  inward  toward  a  more  or  less  uniform  cold  core.  When 
a  sufficient  mass  of  warmed  air  had  accumulated  within  the  thermal  layer,  buoyant  forces 
eventually  broke  up  the  layer  and  allowed  the  unevenly  heated  air  to  stratify  into  hot 
and  cold  layers.  This  idea  was  eventually  rejected  because  the  air  temperature  distributions 
measured  with  the  sliding  probes  showed  no  signs  of  a  thermal  boundary  layer.  A  second 
idea  suggested  that  the  elevated  heat-transfer  rates  might  be  attributable  to  the  roughness 
of  the  inside  of  the  tube.  Accordingly,  the  roughness  of  the  tube  was  measured  and  found 
to  be  about  0.0015  in.  The  resulting  roughness-to-tube  diameter  ratio  indicated  the  tube 
was  essentially  smooth.  Thus,  no  significant  increase  in  heat  transfer  resulted  from  a 
roughness  correction.  It  therefore  appears  that  the  most  feasible  model  of  the  heat-transfer 
process  studied  here  is  that  of  decaying  jet  turbulence  opposing  buoyant  effects  as  discussed 
in  Appendix  B. 
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Several  potentially  useful  relations  have  been  obtained  from  correlating  the  proposed 
model  with  experimental  data.  The  heat-transfer  rate  as  a  function  of  position  in  the 
tube  may  be  computed  from  the  following  empirical  equation,  developed  in  Appendix 
B 

Nud  (x)  =  0.020  [Reo  (1  +  0.0074  e  e-o.0  5  7  T)]  0.87  5 


which  was  determined  from  data  covering  a  diameter  Reynolds  number  range  of  32,000 
to  39,000  and  axial  position  of  0  <  ^  <  88  diameters.  The  Prandtl  number  varied  little 
from  0.7,  and  the  tube-to-jet  area  ratio  (e)  was  353.  The  equation  is  plotted  in  Fig.  12 
with  the  experimental  data  and  the  other  more  conventional  theories  discussed  earher. 
It  is  clear  that  this  relation  correlates  the  data  better  than  the  classical  theories.  As  might 
be  expected  from  the  mechanism  of  decaying  turbulence,  the  proposed  model  is  in  better 
agreement  with  the  free  convection  process. 

An  energy  ratio  correlation  equation  was  developed  in  Appendix  B  (Eq.  (B-22)).  By 
evaluating  the  known  parameter,  Eq.  (B-22))  can  be  rearranged  to  give  the  approximate 
location  of  the  last  unstratified  station. 
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where  the  critical  energy  ratio  was  taken  as  30  from  Fig.  B-2  and  the  Grashof  number 
of  the  data  base  ranged  from  2  x  10^  to  2  x  lO^®.  Uniform  cooling  of  a  tube  of  length 
greater  than  Xcrit  can  be  achieved  by  spacing  multiple  inlet  ports  at  intervals  along  the 
tube  at  spacings  of  about  Xcrit  and  cooUng  the  tube  in  stages.  The  correlation  equation 
also  suggests  the  unstratified  length  of  tube  may  be  extended  by  increasing  the  jet 
turbulence  (reducing  the  jet  flow  area),  though  this  approach  is  practically  limited  by 
the  jet  Mach  number. 


6,0  SUMMARY  OF  RESULTS 

Heat-transfer  rates  measured  in  a  long,  slender,  externally  insulated  pipe  having  forced 
convective  internal  cooling  have  been  found  to  be  much  greater  than  that  expected  from 
classical  theories.  In  addition,  significant  thermal  stratification  has  been  measured  at  various 
distances  from  the  inlet  port,  apparently  depending  upon  the  charge  pressure.  Since  no 
conventional  model  was  found  to  satisfactorily  explain  the  observed  phenomena,  a  model 
was  developed  from  the  basic  idea  of  decaying  turbulence  initially  produced  by  the  inlet 
jet.  In  this  model  the  turbulence  not  only  elevates  the  heat-transfer  rate  near  the  inlet 
but  also  delays  the  onset  of  thermal  stratification.  Decay  of  the  turbulence,  however, 
eventually  allows  free  convection  to  dominate.  The  axial  location  where  stratification  begins 
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(Xcrit)  has  been  correlated  with  a  critical  value  of  the  turbulent-to-buoyant  kinetic  energy 
ratio.  The  gross  effects  of  stratification  and  its  dependence  upon  pressure  have  been  verified 
from  the  present  experimental  results,  but  confirmation  and  refinement  of  the  detailed 
explanation  must  be  based  upon  more  elaborate  testing,  including  measurement  of 
turbulence  along  the  tube  and  measurement  of  the  temperature  distribution  over  the  entire 
cross  section  of  the  tube.  The  results  obtained  so  far,  however,  do  allow  some  estimate 
of  the  effect  of  the  jet  upon  the  heat-transfer  rate  and  upon  the  length  of  tube  that 
may  be  uniformly  cooled  by  depending  upon  jet  turbulence  to  prevent  free  convection. 
In  addition,  it  appears  feasible  to  uniformly  cool  longer  lengths  by  using  multiple  inlet 
ports  to  cool  the  tube  in  stages.  The  distance  between  stages  should  be  about  JCcrit- 
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APPENDIX  A 

COMPUTATION  OF  HEAT-TRANSFER  COEFFICIENTS  FROM 
EXPERIMENTAL  TEMPERATURE  DISTRIBUTIONS 

Two  methods  of  computing  heat-transfer  coefficients  were  used  herein,  one  based 
on  the  time  rate  of  change  of  steel  temperature  at  a  given  axial  location  along  the  tube 
and  the  other  based  on  the  axial  variation  of  the  air  temperature  and  the  mass  flow  rate. 
Consider  the  first  method.  Taking  as  a  definition  for  the  heat-transfer  coefficient  (h)  the 
equation 

q  =  hA0as  (A-1) 

where  q  is  the  heat  transfer  per  unit  surface  area,  and  AO^s  is  the  temperature  difference 
between  the  air  and  steel,  h  can  be  easily  determined  from  the  experimental  data.  For 
a  surface  area  A, 

q  =  Q/A  (A-2) 


The  energy  equation  for  a  control  mass  of  the  steel  yields 

Q  =  Ms  Cp3  [0s(x,t)  -  6ls(x,  t  +  At)]  (A-3) 

where  Ms  is  the  mass  of  the  steel,  Cp^  is  the  specific  heat,  is  the  steel  temperature 
at  a  given  axial  position,  and  t  is  the  time.  For  a  thin  walled  tube, 

Ms  =  Ps  ^tDtAx  (A4) 


where  ps  is  the  steel  density,  D  is  the  diameter  of  the  charge  tube  and  t  its  thickness, 
and  Ax  is  a  short  length  of  tube.  Then,  dividing  Eq.  (A-3)  through  by  At,  noting  that 
A  =  ttDAx,  and  combining  all  these  results  with  the  four  equations,  the  following  data 
reduction  equation  emerges. 


h(x,t*) 


Ps  Cpj  T  [0s(x,t)  ■  0s(x,t  +  At)] 
[0s(x,t*)-0a(x,t*)]  At 


(A-5) 


where  is  the  air  temperature  at  a  given  time  and  position  and  t  <  t*  <  t  +  At.  Recall 
from  the  experimental  data  that  the  steel  temperature  generally  varies  circumferentially 
and  the  air  temperature  vertically.  To  account  for  this,  the  temperatures  in  Eq.  (A-5)  were 
computed  from  weighted  averages  of  all  the  available  data  at  a  given  axial  location.  Both 
the  air  and  steel  temperatures  at  the  extreme  top  and  bottom  of  the  tube  were  weighted 
0.5  while  all  others  were  weighted  1.0.  Another  question  arose  concerning  circumferential 
and  axial  conduction  within  the  tube,  since  these  possibilities  are  not  included  in  the 
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energy  equation  (A-3).  A  manual  calculation  at  the  most  extreme  axial  temperature 
differences  indicated  that  conduction  along  the  tube  was  negligible  compared  to  the  heat 
flux  across  the  film.  On  the  other  hand,  circumferential  conduction  was  found  to  be 
significant,  and  a  precise  accounting  was  worked  out  for  the  energy  equation.  However, 
it  turned  out  that  the  correction  had  no  net  effect  upon  the  weighted  average  steel 
temperature  and  thus  did  not  change  the  experimental  heat-transfer  coefficients  computed 
without  the  correction. 


In  the  second  method,  the  heat-transfer  coefficient  was  computed  on  the  basis  of 
air  temperature  variation.  The  integral  energy  equation  for  a  control  volume  can  be  reduced 
to 


Q  =  m  Cp^  [0a(X2-t)  -  0a(Xl,t)] 


(A-6) 


where  xi  and  X2  are  the  inlet  and  outlet  positions  on  the  control  volume,  Cp^  is  the 
specific  heat  of  air,  and  m  is  the  mass  flow  rate.  Taking  the  area  equal  to  7rD(x2-Xi) 
and  combining  Eqs.  (A-1),  (A-2).  and  (A-6)  produces 


h(x*,t) 


m  Cpa  [^a(x2,t)  ■  0a(xi  ,t)l 

7rD(x2  -  Xi  )  [0s(x*,t)  -  0a(x,*t)] 


(A-7) 


where  xj  <  x*  <  X2.  For  computational  purposes  x*  and  t*  were  taken  as  the  midpoint 
of  the  space  and  time  intervals. 
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APPENDIX  B 

PARTIAL  MODEL  OF  THE  HEAT-TRANSFER  PROCESS 


Several  attempts  were  made  to  model  the  heat-transfer  process  in  terms  of  relations 
which  would  reduce  to  classical  theories.  Three  classical  theories  are  shown  in  Figure  12. 
Relations  are  shown  for  a  turbulent  flat  plate  boundary  layer;  forced  convection  in  fully 
developed,  turbulent  pipe  flow;  and  turbulent  free  convection  inside  a  sphere.  Free 
convection  gives  the  closest  approximation  to  the  experimental  data  but  is  probably  not 
sufficient  to  explain  the  heat-transfer  process  encountered.  Models  in  terms  of  flat  plate 
boundary  layers  and  fully  developed  pipes  modified  by  the  inclusion  of  decaying  jet 
turbulence  were  considered,  but  they  were  both  rejected  because  the  experimental  profiles 
did  not  exhibit  boundary-layer  properties.  The  model  finally  settled  upon  was  derived 
from  a  suggestion  that  the  effect  of  the  jet  turbulence  is  to  reduce  the  boundary-layer 
thickness.  It  was  thus  hypothesized  that  the  jet  turbulence  reduced  the  boundary  layer 
to  essentially  a  laminar  sublayer  so  that  the  heat-transfer  process  became  one  of  pure 
conduction  across  the  layer  with  a  relatively  large  temperature  gradient  varying  from  the 
wall  temperature  to  the  bulk  fluid  temperature  over  the  thickness  of  the  sublayer.  Thus, 
from  the  one-dimensional  Fourier  equation  for  pure  conduction. 


q  =  Ka 


dy 


(B-1) 


where  y  is  the  direction  normal  to  the  tunnel  wall,  and  from  Eq.  (A-1)  the  heat-transfer 
coefficient  may  be  written  as 


Ka 

where  6  is  the  local  sublayer  thickness.  This  result  assumes 

30a  _  AOas 

3y  6 


(B-2) 


(B-3) 


The  following  equation  adequately  describes  the  laminar  sublayer  thickness  for  fully 
developed  pipe  flow 


,  3  :  .8  D 

^  “  Reov/r 

where  f  is  the  friction  factor  given  by 


0.316 

^  ^  2.000 


(B-4) 


(B-5) 
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(see  Ref.  5).  Calculations  based  on  these  equations  give  film  coefficients  consistent  with 
the  experimental  data  not  too  near  the  inlet  port  (see  Fig.  1 2). 


Near  the  inlet  port  where  the  jet  turbulence  is  relatively  intense,  the  sublayer  is 
hypothesized  to  be  of  reduced  thickness  5'  due  to  a  local  turbulence  of  magnitude  u'. 
The  Reynolds  number  in  Eqs.  (B-4)  and  (B-5)  is  thus  written  as 


Reb  = 


p(u  +  u ')  D 
M 


(B-6) 


where  u  is  the  local  bulk  velocity.  Another  more  convenient  form  is 


Combining  Eqs.  (B-4),  (B-5),  and  (B-7)  produces 

6  _  58.3h  D 

^  '  Ka  ReD°-®^  5 


(B-7) 


(B-8) 


Use  of  the  experimental  heat-transfer  coefficients  in  Eq.  (B-8)  then  gives  the  ratio  of 
the  usual  sublayer  thickness  to  that  suggested  by  the  elevated  heat-transfer  rate.  Further, 
from  the  above  equations  is  it  possible  to  write 


u 

u 


(B-9) 


which  gives  an  experimental  estimate  of  the  turbulent  intensity. 


The  above  equations  along  with  classical  Nusselt  number  relations  were  programmed 
for  computer  processing  of  the  experimental  data.  The  statistical  distribution  of  the  u'/u 
data  was  then  examined.  Figure  B-1  shows  the  data  plotted  versus  axial  distance  from 
the  inlet  port.  The  decay  of  the  turbulence  was  correlated  with  distance  from  the  inlet 
with  the  relation 

f 

—  =  2.6  e-o  577  (B-IO) 


where  x  is  the  nondimensional  distance  x/D.  Correlations  with  powers  of  time,  such  as 
for  an  isotropic  turbulence,  were  examined  but  failed  to  correlate  the  data  as  well  as 
Eq.  (B-10).  In  an  effort  to  establish  the  dependence  of  u'/u  on  the  velocity  of  the  inlet 
jet,  the  one-dimensional  continuity  equation  between  the  jet  and  some  uniform  station 
along  the  large  tube  was  invoked. 


PjUjAj  =  puA 


(B-11) 
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For  incompressible  flow 


=  € 


(B-12) 


where  the  j  subscript  refers  to  the  jet  and  e  is  the  area  ratio.  Assuming  then  that 


u 


u 


(B-13) 


at  x/D  =  0  so  that  the  constant  2.6  in  Eq.  (B-10)  actually  should  be  written  as  constant 
times  e,  the  dependence  on  the  diameter  of  the  charging  port  appears  as 

f 

^  =  0.0074  e 

Finally,  defining  a  Nusselt  number  based  on  diameter  as 

/hD\ 

=  fe) 

and  using  Eq.  (B-14)  along  with  the  previous  sublayer  equations,  the  following  Nusselt 
number  relation  is  obtained. 


(B-14) 


(B-15) 


Nud  =  0.020  (Reo  (1  +  0.0074  e  e-0  0  57  x)|0.87  5  prO.4 

To  obtain  this  relation,  Pr®  ^  (for  Pr  =  0.68,  the  experimental  Prandtl  number)  was 
extracted  from  the  numerical  constant  in  order  to  compare  this  equation  with  the  familiar 
relation  for  a  fully  developed  turbulent  pipe  (Ref.  4). 


Nud  =  0.023  Reo^-^  Pr^-^ 

The  multiplicative  constants  are  comparable  in  the  two  relations,  but  the  larger  exponent 
of  Reo  in  Eq.  (B-16)  reflects  the  fact  that  far  downstream  the  heat-transfer  rate  of  this 
experiment  never  decays  back  to  that  of  a  classical  fully  developed  pipe  flow,  the  reason 
being  the  eventual  domination  by  free  convection.  Since  the  asymptotic  value  of  Eq.  (B-16) 
for  large  x  should,  therefore,  be  a  function  of  the  strength  of  the  free  convection,  the 
relation  must  be  considered  deficient  for  its  lack  of  dependence  upon  Grashof  number. 
This  obvious  extension  was  not  attempted  because  of  the  narrow  range  of  Grashof  number 
covered  by  the  experimental  data.  The  relation  given  by  Eq.  (B-16)  is  shown  in  Fig.  12 
with  the  experimental  data  from  which  it  was  derived. 

Before  proceeding,  it  should  be  understood  that  the  replacement  of  Rcq  by  an 
effective  turbulent  Reynolds  number  Reb  in  Eq.  (B-4)  to  obtain  a  relation  for  6',  the 
effective  sublayer  thickness  in  the  presence  of  turbulence,  is  not  rigorous  because  such 
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an  application  of  Eq.  (B-4)  is  probably  not  consistent  with  the  logic  behind  its  classical 
development.  However,  note  from  Eq.  (B-4)  that  increasing  the  Reynolds  number  decreases 
the  thickness  of  the  laminar  sublayer,  a  result  which  would  increase  the  heat-transfer  rate 
for  a  given  temperature  difference  across  the  layer.  Since  this  agrees  qualitatively  with 
the  experiment,  it  is  suggested  that  the  jet  turbulence  increases  the  effective  Reynolds 
number.  A  second  point  to  consider  is  that  in  Eq,  (B-9)  the  relative  turbulence-induced 
change  in  the  layer  thickness  depends  not  on  u'  but  on  the  ratio  of  the  turbulent  intensity 
to  the  local  bulk  velocity  through  the  tube.  For  the  nominally  constant  mass  flow  rates 

of  the  experiments,  the  jet-to-bulk  velocity  ratio  is  also  nominally  constant,  according 

to  the  assumptions  inherent  in  Eq.  (B-11).  The  logical  conclusion,  given  in  these  two 

points,  is  that  the  heat-transfer  rate  should  not  vary  significantly  from  experiment  to 

experiment  (varying  only  pressure  and  not  mass  flow  rate).  This  is  consistent  with  the 
experiment  because  within  the  scatter  bands  of  Fig,  12,  no  consistent  variation  with 
pressure  was  observed  (changing  the  pressure  changes  u  and  uj  but  not  uj/u  according 
to  continuity).  Thus  the  expedient,  but  nonrigorous,  substitution  is  justified  because  it 
is  qualitatively  consistent  with  experimental  observation.  Obviously,  the  model  to  this  point 
would  be  much  stronger  if  an  experimental  confirmation  of  the  dependence  of  the  laminar 
sublayer  thickness  on  free-stream  turbulence  were  available. 

Besides  the  heat-transfer  model  just  described,  further  data  analysis  was  performed 
in  an  effort  to  understand  the  stratification  in  the  air  temperature,  which  was  manifested 
in  terms  of  circumferential  gradients  in  the  steel  temperature.  The  resulting  model  has 
evolved  into  the  following  form. 

From  the  experimental  data  the  following  facts  have  been  established. 

1 .  Circumferential  temperature  gradients  appear  in  the  steel  only  when  the 
air  has  stratified  into  hot  and  cold  layers  at  the  top  and  bottom  of  the 
charge  tube. 

2.  The  onset  as  well  as  the  rate  of  increase  of  stratification  is  delayed  more 
and  more  as  the  pressure  is  lowered. 

Stratification  -  vertical  temperature  gradients  -  is  associated  with  density  gradients;  hot 
air  rises  while  coot  air  falls.  But  stratification  does  not  occur  as  long  as  there  is  sufficient 
turbulence  to  prevent  it.  It  is  therefore  suggested  that  the  onset  of  stratification  can  be 
correlated  with  an  energy  ratio  (er)  equal  to  the  turbulent  kinetic  energy  divided  by  the 
buoyant  kinetic  energy. 


er 


(B-17) 
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When  the  jet  turbulence  has  decayed  sufficiently,  a  small,  warm  lump  of  fluid  with  velocity 
Ub  manages  to  retain  its  identity  long  enough  to  migrate  to  the  top  of  the  tube  where 
its  higher  temperature  reduces  the  heat-transfer  rate  from  that  occurring  at  the  bottom 
of  the  tube.  When  a  sufficient  number  of  hot  lumps  survive  the  migration  time,  gradients 
in  the  steel  begin  to  appear.  The  onset  of  stratification  is  thus  to  be  correlated  with 
the  energy  ratio  becoming  small  enough,  say  less  than  ercrit- 

To  compute  er  from  experimental  data,  Eq.  (B-14)  was  used  for  the  numerator  of 
Eq.  (B-17).  The  buoyant  velocity  Ub  was  computed  by  equating  the  buoyant  force  on 
an  element  of  fluid  with  the  drag  produced  at  the  terminal  velocity  allowed  by  viscous 
drag.  This  approach  was  taken  because  the  unrestrained  acceleration  attributable  to  the 
buoyant  force  produced  absurdly  large  velocities  and  short  migration  times.  From 
Archimedes'  principle  the  buoyant  force  less  the  weight  of  the  element  may  be  written 
as 


Fb  =  (p  -  p')g  V  (B-18) 

where  V  is  the  volume  of  the  element,  p  is  its  density,  g  is  the  acceleration  of  gravity, 
and  p  is  the  density  of  the  surrounding  fluid.  The  drag  on  the  element  can  be  written 
as 

Fd  =  T  C'd  P  Ub"  A  {B-19) 

where  Cp  is  a  drag  coefficient  and  A  is  the  projected  area  of  the  element.  The  use  of 
Eq.  {B-19),  which  represents  drag  on  a  rigid  body,  requires  further  discussion.  An  alternate 
assumption  would  be  to  represent  the  drag  in  terms  of  a  viscous  shear  layer  developing 
between  the  element  and  the  surrounding  fluid  (a  "nondiscrete"  element  approach).  Such 
an  assumption  could  take  the  form  of 

Fd  =  P  Ub  A/nd  (B-19a) 

where  nd  represents  the  thickness  of  the  shear  layer.  Underlying  the  choice  between  Eq. 
(B-19)  and  (B-19a)  is  the  question  of  whether  the  elements  really  are  discrete.  If  Eq. 
(B-19a)  is  the  appropriate  model,  then  the  heated  element  must  be  losing  significant  mass 
to  its  surroundings  via  the  shear  layer  and  is  not  discrete.  If,  however,  Eq.  (B-19)  is  more 
appropriate,  then  the  shear  layer  is  very  steep  and  approximates  a  more  or  less  solid 
boundary,  a  discrete  mass.  Obviously,  neither  equation  can  rigorously  represent  the  entire 
model.  In  defense  of  the  discrete  model,  however,  it  can  be  said  that  the  driving  force 
(buoyancy,  which  is  really  a  body  force)  acts  uniformly  on  all  portions  of  the  element 
to  make  them  move  together.  But  the  basic  logic  of  the  discrete  assumption  is  that  if 
the  heated  elements  do  not  remain  sufficiently  discrete  for  a  significant  duration  of  time 
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SO  as  to  transport  vertically  a  significant  mass  of  heated  air,  then  classical  free  convection 
would  not  occur.  Proceeding  with  the  discrete  element  theory,  some  decision  on  geometry 
of  the  element  must  be  made.  The  approach  taken  here  has  been  to  idealize  the  element 
as  a  discrete  piece  of  fluid  of  constant  geometry  -  namely,  a  sphere  of  diameter  d  -  for 
which  A  =  7rd^  and  V  =  l/brrd^.  With  these  relations  and  the  ideal  gas  law  used  in  the 
relation  resulting  from  equating  Fb  and  Fq  ,  the  following  equation  emerges. 


Ub 

u 


0.31 


GrpO.564 

Reo 


(B-20) 


where  Cq  =  2.67  was  fitted  to  the  drag  coefficient  for  a  sphere  (Ref.  5). 

Here,  Grp  is  the  Grashof  number  defined  as 


Gro  = 


p2  g  A0 
^20 


(B-21) 


and  d/D  is  the  ratio  of  the  element  size  to  the  tube  size,  clearly  a  variable  quantity  since 
any  hot  element  formed  would  dissipate  heat  and  be  gradually  reduced  in  size  by  shear 
stresses.  However,  to  adopt  some  hopefully  typical  value,  d/D  was  taken  to  be  of  the 
order  of  6  -  6',  the  portion  of  the  laminar  sublayer  torn  away  by  the  turbulence.  This 
calculation  gave  an  estimate  of  d/D  =  0.01,  which  produced  reasonable  terminal  velocities 
Ub  and  corresponding  migration  times.  Finally,  substitution  of  Eqs.  (B-14)  and  (B-20) 
in  the  definition  of  the  energy  ratio  Eq.  (B-17)  yields 


er 


0.024  e  Reo  e 


(B-22) 


Figure  B-2  shows  the  energy  ratio  plotted  versus  axial  station.  The  computation  was  made 
using  Rep  and  Grp  computed  from  experimental  data.  The  bands  are  spaced  lone  standard 
deviation  and  represent  scatter  in  the  temperature  distribution  as  well  as  its  time-wise 
variation.  The  parallelograms  locate  the  last  unstratified  stations  within  the  distance  of 
two  succeeding  measurement  locations.  No  location  is  shown  for  the  400-psig  case  because 
the  flow  was  already  heavily  stratified  as  the  first  measurement  location  at  7.5  diameters. 
The  value  of  erciit  in  Fig.  B-2  would  appear  to  be  nominally  of  order  lO^.  The  model 
thus  suggests  that  stratification  begins  when  the  energy  ratio  drops  below  this  value. 


An  approximate  test  of  Eq.  (B-22)  may  be  performed  by  noting  that 

er  a  p-i/4  e-0-ii4  3r  (B-23) 

since  Rep  «  p  and  Grp  «  p2  while  p  a  P.  Taking  the  natural  log  of  Eq.  (B-23)  and 
evaluating  at  two  sets  of  conditions  for  a  constant  energy  ratio  produces 
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X2  -  xj  =  2.2  Cn  — —  (B-24) 

92 

which  is  the  difference  in  axial  location  at  which  stratification  should  occur.  Comparing 
the  1  50-psig  case  with  75  psig  yields  a  difference  of  about  two  diameters  while  comparing 
Figs.  7  and  8  yields  four  diameters,  nominally  the  same  result  considering  the  crudeness 
of  the  model. 

It  is  important  to  note  at  this  point,  should  the  fact  have  been  obscured  by  all 
the  detail,  that  neither  the  agreement  between  the  model  and  the  experimental  stratification 
point  nor  the  agreement  in  Fig.  12  between  the  heat-transfer  model  and  the  data  provides 
a  sufficient  verification  of  the  proposed  model.  The  comparison  does  illustrate  that  the 
model  exhibits  the  gross  level  trends  of  the  data,  but  the  degree  of  agreement  is  indicative 
only  of  the  correctness  of  the  calculations  to  determine  the  model's  constants  from  the 
data. 

In  summary,  a  model  of  the  heat-transfer  process  and  stratification  development  has 
been  established  from  the  idea  of  jet  turbulence  opposing  buoyancy.  The  obvious  weak 
points  in  the  model  are  the  reduced  sublayer  and  the  discrete  element  drag  assumptions, 
as  both  are  based  on  essentially  misused  empirical  relations.  However,  in  defense  of  the 
present  modeling  effort,  it  can  be  said  that  the  original  goal  of  the  effort  was  to  obtain 
some  sensible  governing  equation  by  carrying  through  from  first  principles  the  gross  level 
effects  postulated  to  be  governing  the  present  physical  situation.  Obviously,  improvement 
of  the  turbulence-buoyancy  model  or  derivation  of  a  better  one  would  require  a 
considerable  amount  of  further  experimental  work. 


Distribution  Peaks  and 
Scatter  Bands 


Axial  Station.  x/D 

Figure  B-1.  Decay  of  jet  turbulence. 
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Figure  B-2.  Energy  ratio  versus  axial  station. 
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NOMENCLATURE 

Area,  cross  sectional  or  surface  area  depending  upon  use 
Drag  coefficient 

Specific  heat  at  constant  pressure 
Tube  inside  diameter 
Element  diameter 

Ratio  of  turbulent  to  buoyant  kinetic  energy 
Buoyant  force 
Drag  force 
Friction  factor 

Grashof  number  based  on  tube  diameter 

Acceleration  of  gravity 

Heat-transfer  coefficient 

Thermal  conductivity 

Mass 

Mass  fiow  rate 

Nusselt  number  based  on  diameter 

Number  of  elemental  diameters  in  the  sliear  layer 

Press  u  re 

Prandtl  number,  ^Cp/K 

Net  heat  transferred  (not  heat-transfer  rate) 

Heat-transfer  rate 
Heat-transfer  rate  per  unit  area 
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Rep  Reynolds  number  based  on  tube  diameter 

Rej  Reynolds  number  based  on  size  of  the  heated  element  of  tluid 

t  Time 

Li  Velocity 

V  Volume 

X  Axial  distance  along  tube 

y  Normal  distance  from  tube  wall 

5  Sublayer  thickness 

e  Area  ratio.  A/Aj 

6  Temperature 

p  Mass  density 

p  Dynamic  viscosity 

r  Tube  wall  thickness 

SUBSCRIPTS 
a  Air 

b  Buoyant 

crit  Critical  value  at  which  stratification  begins 

j  Jet 

s  Steel  (or  tube) 

as  Air-steel,  as  in  A0as.  the  air-steel  temperature  difference 

1,2  Initial  and  final  stations  of  a  control  volume  taken  from  the  tube 
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SPECIAL  SYMBOLS 

~  (over  tilda)  Nondimensionalization.  as  in  x  =  x/D 

(over  dot)  Time  derivative,  as  in  Q,  the  net  heat-transfer  rate 
*  Intermediate  value,  as  in  t  <  t*  <  t  +  At 

A  Increment,  as  in  At 

'  (prime)  Turbulent  quantity,  as  in  u',  or  quantity  modified  by  the  presence 

of  turbulence,  as  in  6',  Ren,  or  p' 

oc  Proportional  to 
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